WHAT IS CLAIMED IS: 



1 . A method of modeling the behavior of a molecule, comprising 
selecting a torsion angle, rigid multibody model for said molecule, said model 

having equations of motion; 

selecting an implicit integrator; and 

generating an analytic Jacobian for said implicit integrator to integrate said 
equations of motion so as to obtain calculations of said behavior of said molecule. 

2. The method of claim 1 wherein said analytic Jacobian is derived from 
an analytic Jacobian of the Residual Form of the equations of motion. 

3. The method of claim 2 wherein said anal3^ic Jacobian J comprises 
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where q are the generalized coordinates, u are the generalized speeds. Wis a joint map matrix 
and Mis the mass matrix and p„ is the dynamic residual of the equations of motion, and z is 
-M->„(g,M,0). 

4. The method of claim 3 wherein said impUcit integrator selecting step 
comprises an L-stable integrator. 



5. A method of simulating the behavior of a physical system, comprising 
modeling said physical system with a torsion angle, rigid multibody model, 
said model having equations of motion; and 
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integrating said equations of motion with an implicit integrator; said implicit 
integrator having an analytic Jacobian to obtain calculations of said behavior of said physical 
system. 

6. The method of claim 5 wherein said analytic Jacobian is derived from 
an analytic Jacobian of the Residual Form of the equations of motion. 

7. The method of claim 6 wherein said analytic Jacobian J comprises 
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where q are the generalized coordinates, u are the generalized speed, Wis a. 
joint map matrix and Mis the mass matrix and /)„ is the dynamic residual of the equations of 
motion, and z is -M~^p^(q,u,0). 

8 . The method of claim 7 wherein said implicit integrator comprises an 
L-stable integrator. 

9. Computer code for simulating the behavior of a molecule, said code 

comprising 

a first module for a torsion angle, rigid multibody model of said molecule, said 
model having equations of motion; and 

a second module for an implicit integrator to integrate said equations of 
motion with an analytic Jacobian to obtain calculations of said behavior of said molecule. 

10. The computer code of claim 9 wherein said analytic Jacobian is 
derived from an analj^ic Jacobian of the Residual Form of the equations of motion. 
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1 1 . The computer code of claim 1 0 wherein said analytic Jacobian J 



comprises 
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where q are the generalized coordinates, u are the generahzed speed, W'l&z 
joint map matrix and M is the mass matrix and yO„ is the dynamic residual of the equations of 
motion, and z is -M'^p^ (q,u, 0) . 

12. The computer code of claim 1 1 wherein said implicit integrator 
comprises an L-stable integrator. 

1 3 . Computer code for simulating the behavior of a physical system, said 
code comprising 

a first module for a torsion angle, rigid multibody model of said system, said 
model having equations of motion; and 

a second module for an implicit integrator to integrate said equations of 
motion with an analytic Jacobian to obtain calculations of said behavior of said system. 

14. The computer code of claim 13 wherein said analj^c Jacobian is 
derived from an analytic Jacobimi of the Residual Form of the equations of motion, 

15. The computer code of claim 14 wherein said analytic Jacobian J 

comprises 
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where q are the generahzed coordinates, u are the generalized speed, is a 
joint map matrix and M is the mass matrix and p„ is the dynamic residual of the equations of 
motion, and z is -M~^p^(q,u,0) . 

16. The computer code of claim 1 5 wherein said implicit integrator 
comprises an L-stable integrator. 
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